Josephson oscillations in binary mixtures of F = 1 spinor BECs 
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We analyze theoretically Josephson oscillations in a mixture of two Zeeman states of a spinor 
Bose-Einstein condensate in a double-well potential. We find that in the strongly polarized case, 
the less populated component exhibits a complex dynamics with an anti- Josephson behavior, i.e. 
oscillates in phase with the more populated one. In the balanced population case the Josephson 
oscillations unveal a dependence with the different spin collision channels. This effect could be used 
to experimentally measure the distinct scattering lengths entering in the description of a spinor 
condensate. Our numerical results are in close agreement with an analytical description of the 
binary mixture using a two mode model. 
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Tunneling is one of the most fascinating quantum phe- 
nomena known since the development of quantum me- 
chanics. At macroscopic scale the Josephson effect essen- 
tially consists on the fast oscillating tunneling through a 
macroscopic barrier driven by a quantum phase difference 
between the two sides of the potential barrier [ij. 

Josephson dynamics of individual atoms has been the- 
oretically studied by several groups Its extension to 
the tunneling of a Bose-Einstein Condensate (BEC) in a 
two- well potential was presented in On the experi- 
mental side, Josephson tunneling of BEC in an optical 
lattice was first reported in Ref. Recently, a clear 
evidence of a bosonic Josephson junction in a weakly 
linked scalar BEC has been presented [H. There, the 
observed time evolution of both, the phase difference be- 
tween the two sides of the barrier and the population 
imbalance shows the importance of the interactions and 
is accurately described by directly solving the full Gross- 
Pitaevskii equation (CP) or by using a two-mode approx- 
imation 0. 

Multi-component BECs trapped in double wells of- 
fer an interesting extension of the tunneling problem. 
Binary mixtures (i.e. pseudo-spin 1/2 BECs) that sup- 
port density-density interactions have been analyzed in 
Refs. 0, using mainly the two-mode approach. More 
complex spinor BECs, that support population transfer 
between different Zeeman states have been also studied 
as possible atomic candidates for the macroscopic quan- 
tum tunneling of mag netization (MQTM) [i,[i3|. 

In this letter, we investigate the dynamics of a weakly 
linked mixture of m — ±1 Zeeman components of a = 1 
spinor ®'^Rb condensate. By solving the CP equation, we 
find that in the strongly polarized case the less populated 
component oscillates in phase with the more populated 
one. For equal populations, the corresponding Joseph- 
son frequency provides information about the different 



scattering lengths of the system. We show how these 
two combined setups give access to the spin dependent 
interaction strenght. Our method extracts information 
from the macroscopic dynamics of the two Zeeman com- 
ponents on a two-well potential. It can be seen as the 
macroscopic counterpart of the method of Ref. [11] where 
the spin dynamics of many atom pairs is considered based 
on an earlier observation of coherent spin-dynamics of ul- 
tracold atom pairs trapped in sites of an optical lattice in 
a Mott insulator regime. Our method would improve the 
current knowledge, which despite the accuracy of their 
involved experiments cannot, for instance, uniquely de- 
termine whether the F — 2 ground state of ^^Rb is in an 
anti- ferromagnetic or a cyclic phase 11 1 . 

Following the conditions used in the experiments of 
Josephson tunneling in a scalar condensate we con- 
sider N — 1150 atoms of spin-1 ^^Rb trapped initially in 
the following potential: 

F(r) = {M/2){colx^+u;y+u;y)+VoCOs''{7:{x~Am)/qo), 



with ( 



2tt X 78 Hz, w„ = 27r X 66 Hz, = 27r x 90 Hz, 



qo = 5.2 /Ltm, Vb — 413 /i Hz, and M is the mass of the 
atom. Am defines the initial asymmetry in the double- 
well potential and thus is directly related to the initial 
populations imbalance. At t = the small asymmetry 
is switched off (A^ = 0) and the system is let to evolve 
in the double-well potential, which is the same for all 
components. The fact that the double well is created 
only along the x-direction forces the dynamics of the 
Josephson oscillation to be one dimensional as beautifully 
demonstrated in Then, defining u!± — ^JuJ^uJ^ the 
couphng constants can be rescaled by a factor l/(27ra^), 
with a_i_ the transverse oscillator length and the dynam- 
ical equations transform into one-dimensional ones d2j|. 

In the mean-field approach, the F = 1 spinor conden- 
sate is described by a vector order parameter ^' whose 
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components V'm = |'0m|e*'''"' correspond to the wave func- 
tion of each magnetic sublevel \F = l,m) = |m) with 
TO = 1,0, — 1. In absence of an external magnetic field 
and at zero temperature the spin dynamics of this system 
confined in an external potential, V, is described by the 
following coupled equations for the spin components [l3| : 

= [Ws + C2(n±i + no - n=pi)]V'±i + C2V'oV'j^i ; 
^ [Us + C2{ni + n^i)]ipo + C22ipiipQ->p_i , (1) 

with Hs = -h^/{2M) V'^ + V + cqu being the spin- 
independent part of the Hamiltonian. The density of 
the TO-th component is given by ^^(r) = j^/jmCr)!^, while 
^(i") = J2m is the total density normalized to 

the total number of atoms N. The population of each 
hyperfine state is Nm = J dr|i/;m(r)p. The couphngs are 
Co = 477^2 (ao-h2a2)/(3M) and C2 = 47rfi,2(a2-ao)/(3M), 
where oo and 02 are the scattering lengths describing 
binary elastic collisions in the channels of total spin 
and 2, respectively. Their values are ao — 101. Sas and 
02 = 100.4aB[l3]. The interatomic interactions permit 
the transfer of population between the different Zeeman 
components by processes that conserve the total spin, 
|0) -I- |0) ^ I + 1) + I — 1). Therefore, in a spinor conden- 
sate the number of atoms in each component is not a con- 
served quantity. Due to the chosen conditions, notably 
the small total number of atoms and the fact that we will 
consider an initial zero population of m = atoms, the 
population transfer effects between the different compo- 
nents is very small 1^. Therefore, in our calculation the 
number of atoms in each Zeeman sublevel remains con- 
stant in time as in a real binary mixture, but the spinor 
character is preserved through the parameter C2 which is 
the spin-dependent collision term. 

In the considered double-well potential we can de- 
fine the population imbalances between the two sides of 
the trap for each component as: Zm{t) = [Nm,L{'t) — 
Nm,R{'t)]/Nm{t) where N^i^(^fi'^{t) corresponds to the 
population on the left (right) side of the trap. The other 
relevant quantity is the phase difference between both 
sides of the trap defined as, (pm{t) — (pm,Rit) — <^m,L(i) 
for each component. The GP equation predicts the evo- 
lution of function of both t and x. Within our 
parameters, the condensates at each side of the trap re- 
main mostly in a coherent state meaning that the x de- 
pendence of (fm is fairly small, as shown in Fig. [T] where 
the condensates in each potential well have an almost 
constant phase. Thus, for simplicity, we define 'fm,L(R) 
as the phase of ipmixM), where xm is the position of the 
maximum density of each condensate. 

For the sake of clarity we start by analyzing as a ref- 
erence the scalar Josephson effect. To this aim, we put 
all the atoms in a single component and solve the corre- 
sponding time dependent GP equation equivalent to the 
scalar case with a non linear term proportional to C0-I-C2. 




FIG. 1: Time evolution of the population imbalance (dashed) 
and phase difference (solid) for a scalar BEC. The upper left 
inset depicts a snapshot of the population density at t=0 ms 
(solid) and t = 10 ms (dashed). The upper right inset shows 
the phase as a function oi x a,t t — 120 ms. The solid line 
stands for the scalar BEC with initial population imbalance 
0.06. The dashed and dotted lines correspond, respectively, 
to the m = +1 and m = —1 components of a simulation with 
{Ni/N,No/N,N-i/N) = (5%,0%,95%) and 2±i(0) = tO.06. 
In all cases the initial phase is set to zero for all x. 



In Fig. [T] we show the time evolution of the population 
imbalance and phase difference for A — — 0.20/im. Both 
quantities oscillate with a frequency luj which slightly dif- 
fers from the bare Rabi frequency (noninteracting limit) 
by a term proportional to the interactions [sj . Two snap- 
shots of the density profile at i = and i = 10 ms can 
be seen in the inset. 

New interesting phenomena appear when a binary mix- 
ture is considered. The binary mixture is simulated by 
means of Eqs. ([T]) with an initial zero amount of to = 
atoms. We have checked however, that all the results pre- 
sented in this paper are nevertheless stable against the 
presence of a small amount of m = atoms in the trap 
{Nq/N < 1%). We consider A±i = =F0.10^m, which pro- 
duces asymmetric initial density profiles for both compo- 
nents with imbalances of 2:±i(0) = ^0.06. The configu- 
ration considered is N-i/N = 95% and Ni/N = 5% and 
no initial phase difference. 

The system evolves in the following way. The most 
populated component essentially drives the dynamics, its 
population imbalance and phase difference, see Fig. [21 fol- 
low a behavior qualitatively and quantitatively similar to 
that of the scalar condensate. Whereas the less populated 
component, to = -1-1 in this case, follows the dynamics of 
the other one. First it does an "anti- Josephson" tunnel- 
ing, i.e. the absolute value of its population imbalance 
actually grows driven by the other component, contrary 
to what would have happened in the absence of the other 
component. This is essentially due to the prevalence of 
the interaction with the other component over the popu- 
lation imbalance driven Josephson tunneling. In addition 
to the high frequency Josephson oscillation, tuj, the less 
populated component has an additional lower frequency 
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total number of atoms of each species, can be written as: 
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FIG. 2: Population imbalance (upper panel) and phase differ- 
ence (lower panel) between the two sides of each condensate. 
The dotted (black) line corresponds to the scalar case, the 
solid and dashed blue lines correspond to the m = +1 and 
m = — 1 components of a simulation with (5%, 0%, 95%). The 
red dot-dashed lines correspond to the (50%, 0%, 50%) case. 
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FIG. 3: Population imbalance (upper panel) and phase dif- 
ference (lower panel) between the two sides of each BEC for a 
simulation with (50%, 0%, 50%). Solid (red) lines correspond 
to the m = — 1 atoms and dashed (blue) to the m=-|-l atoms. 



oscillation with a period of nearly 150 ms. Both effects 
will be analyzed below within the two mode approxima- 
tion. 

A second interesting configuration is that of having 
the same population for both m = ±1 components, 
(50%,0%,50%) with the same initial population imbal- 
ance as before, z±i{0) = =f0.06. In this case the behav- 
ior is completely different. Essentially, as can be seen in 
Fig. [21 a Josephson tunneling with an oscillation period 
of approximately 150 ms (almost 6 times larger than in 
the scalar case) is observed for both components. Thus, 
the longer oscillation, which actually corresponds to the 
lower frequency oscillation seen in Fig. [2l is a pure effect 
of having a binary mixture. 

To get a deeper insight on the observed behavior we 
have performed a two mode analysis of a binary mixture 
of BECs. The general coupled GP equations for a binary 
mixture of components a and b, which keeps invariant the 
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These general couplings, gij, can be related to the cq and 
C2 ones entering in Eqs. ([T]) as gaa — gbb = cq + C2 and 
gab — gba = Cq — C2- Since for *^Rb atoms in F = 1 
|c2| << Co, the mixture we are considering has the same 
self interaction between each component and similar to 
the mixing term, gab ~ gaa- 

The two mode approximation can be expressed with 
the following ansatz 3]: ^'^ = jL{t)'^jL{r) + 
-^ jR{t)<^jR {r), with {<^^a\'^jp) = (5y5a/3 and ^'^■^(i) = 
\/iV,-,a(i)e"^^'°(*', i,j = a,b and a,/3 = L,R. Equations 
for the coefficients ^^.^(t) are obtained from Eqs. ^ ne- 
glecting terms involving mixed products of L and R wave 
functions of order larger than 1. Rewritten in terms of the 
imbalances Zj{t) — [Nj^i^{t) — Nj^^{t)]/Nj{t), and phase 



differences, 



t>j^L they take the form 
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where uj. 
K = - I df 

Uij = gij 



VI -b 

2K/h denotes the Rabi frequency, and 



(3) 



v$L(fO • v$fl(f) + $i(f) i/$fl,(f) 
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Lir)<i>lLir) ^ g^J 



dr'^>lR{r)<^lR{r) . 



In Eqs. ([3]) we have considered a symmetric double- well 
and the same mass for the particles of both components 
TOa = TOfc = M, thus defining ^l{r) = ^a,L{R) = ^b,L{R)- 
This yields: Uaa = Ubb = U and Uab = Uba = U. The 
stability of these equations has been analyzed in Ref . @] . 

The dynamics obtained by solving the coupled GP 
Eqs. Il]) can be now understood by taking the appro- 
priate limits on Eqs. ([3]), after linearizing them around 
the equilibrium values Za = Zb = = (pa — (f'b- Defining 
a and f3 as Na = N{1 - a), Nb = Na and U = U{1 + /3), 
Eqs. dSI reduce to 

Za = ~i-^r4'a , 4 = -t^r06 

4 = [{l-a)NU/h + LUr]za + a{l+ f3)NU/hzb 
4>b = il-a){l + f3)NU/hza + [aNU/h + iOr]zb.i4:) 

These equations correspond to two coupled nonrigid pen- 
dulums with normal mode frequencies lu^ and w<. In 
our system, gab ~ gaa — gbb, i-e. U ^ U and thus 
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f3 « 1. 

by: w> = 



Then the eigenfrequencies can be approximated 

The behavior observed in Fig. [5] is explained by taking 
a — > 0. In this hmit the frequencies simphfy to a;> = LUr 
and cj< = ujr + UN/{hujr) = Loj . The highest popu- 
lated component, Za, decouples from the less populated 
and thus performs a Josephson oscillation, as shown in 
Fig. [51 with the same frequency as the scalar case uoj. 
And the less populated component oscillates with the 
two eigenfrequencies: ujy — Ur and uj^ — ojj. 

The case of equal populations can also be addressed 
within Eqs. ([4|) by setting a = 1/2. In this case the long 
mode can be enhanced by starting with opposite initial 
imbalances as in Fig. [31 Then, both components oscillate 



with 



NU 



(3 = Co. Since /3 << 1, this explains 



the similarity between the oscillation frequency of the 
(50%,0%,50%) case studied and the long oscillation of 
the less populated component of the (5%, 0%, 95%) case, 
see Fig. [51 Let us emphasize that the frequency shift 
with respect to the Rabi frequency, NU / {4:h)i3, is a direct 
consequence of the difference between gaa and gat, which 
in turn depend on the interaction strengths, cq and C2. 
Alternatively, if the scattering lengths are known, the 
Josephson oscillation provides a tool to determine with 
high precision the number of atoms. 

It is important to mention that we have assumed zero 
temperature and therefore, no phase fluctuations are 
present in our calculation. Indeed, phase fluctuations will 
evidently destroy the phase coherence of the condensate 
in each potential well leading to a complete suppression 
of Josephson oscillations. For very elongated systems, it 
is known that even at small temperatures the phase of 
the condensate fluctuates [3]. However, in the setup we 
have considered, which is closely related to the experi- 
mental setup demonstrating Josephson oscillations in a 
scalar BEC 's'l , phase fluctuations due to thermal effects 
will remain small since the condensate is essentially three 
dimensional although the dynamics is one dimensional. 
A detailed analysis of possible noise sources using a two 
mode approximation has been developed for a scalar con- 
densate in a double well potential in showing that 
coherence is preserved at the temperatures at which the 
experiment runs, and demonstrating also how coherence 
is lost by increasing the temperature. A similar analysis 
applies here with a spinor BEC. 

Summarizing, we have addressed the problem of a bi- 
nary mixture of BECs which is realized with a _F = 1 
spinor BEC where the m — component is not popu- 
lated. The behavior of the mixture has been analyzed in 
two significant limits: a binary mixture with a large dif- 
ference in the concentration of the two components and a 
binary mixture with equal amount of both components. 
Both cases have been studied with the CP equations and 
have also been scrutinized by employing a simple two 
mode approximation. 



An independent measurement of the Rabi frequency 
could be obtained from the long oscillation of the less 
populated component of Fig.[2l That would allow for pre- 
cise determination of the (otherwise difficult to measure, 
cf. I 111] ) spin-dependent interaction strength, 02 — 09, us- 
ing the relations, (02 — ao)/{ao + 202) — — /3/(2 + (i) and 
[3 can be evaluated by using the above frequencies as: 
P = iujr {uJr — / {ojj — (jJr)- The accuracy of the method 
is expected to be comparable to that of Ref. pj| and is 
only constrained by the ability to measure with precision 
the different macroscopic oscillations. The method pre- 
sented could be extended to any atom species and would 
be easy to generalize to the F = 2 case. 

We thank M. Oberthaler for encouraging us to develop 
this research direction. B.J-D. is supported by a CPAN 
CSD 2007-0042 contract. This work is also supported 
by Grants No. FIS2008-01661, FIS2008-00421, FIS2008- 
00784, FerMix, and QOIT from MFC. 

Note added. After the submission of this paper, a re- 
lated study appeared, reporting a two-mode study of bi- 
nary mixtures [l6| . 
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